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ABSTRACT 

We investigate the consequences of a continuously injecting central engine on the gamma-ray burst 
afterglow emission, focusing more specifically on a highly-magnetized millisecond pulsar engine. For 
initial pulsar parameters within a certain region of the parameter space, the afterglow lightcurves are 
predicted to show a distinctive achromatic bump feature, the onset and duration of which range from 
minutes to months, depending on the pulsar and the fireball parameters. The detection of or upper limits 
on such features would provide constraints on the burst progenitor and on magnetar-like central engine 
models. An achromatic bump such as that in GRB 000301C afterglow may be caused by a millisecond 
pulsar with P = 3.4ms and B p = 2.7 x 10 14 G. 
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1. INTRODUCTION 

Much of the current research on Gamma-Ray Bursts 
(GRB) is aimed at determining the nature of the central 
engine and its progenitor system. While recently substan- 
tial results have begun to accumulate, the evidence is still 
tentative. Thus, the investigation of criteria to differenti- 
ate between various central engine possibilities is desirable. 
Almost all the present fireball models, including those con- 
sidering various non-uniform injection scenarios, assume 
that the energy injection into the fireball occurs in a short 
period of time. This is also the case in the "refreshed 
shock" scenario (Rees & Meszaros 1998; Kumar & Piran 
2000; Sari & Meszaros 2000). However, in some types 
of central engines, such as a fast-rotating high-field pulsar 
(magnetar) or a black hole plus a long-lived debris torus 
system, a significant energy input into the fireball may in 
principle continue for a time scale significantly longer than 
the 7-ray emission. Therefore, there is a need to investi- 
gate a continuously-fed fireball in more detail. An addi- 
tional motivation is provided by the recent detection of Fc 
features in the X-ray afterglow of GRB 991216 after about 
1.5 days (Piro et al. 2000) and GRB 000214 after about 
1 day (Antonclli et al, 2000), which may require a contin- 
uing post-burst outflow in order to achieve less restrictive 
Fe abundance constraints (Rees & Meszaros, 2000). Dai 
& Lu (1998a, b) first considered continuous injection from 
a millisecond pulsar to interpret the afterglow light curves 
of some GRBs, but did not perform a systematic study 
on this topic. In this paper, we investigate the observa- 
tional consequences of a continuously injecting central en- 
gine, and more specifically, focus on the possibility that 
the central engine is a millisecond pulsar, in particular a 
highly magnetized pulsar or magnetar. 

2. CONTINUOUS-INJECTION DYNAMICS 

We consider an engine which emits both an initial impul- 
sive energy input E- lmp as well as a continuous luminosity, 
the latter varying as a power law in the emission time. In 
this case a self-similar blast wave is expected to form at 



late times (Blandford & McKee, 1976). The differential 
energy conservation relation for the self-similar blast wave 
can be written as dE/dt = C (t/t ) q — n'(E/t), where E 
and t are the energy and time measured in the fixed frame. 
The first term £ = Co(t/to) q , where q' and k' are dimen- 
sionless constants, denotes the continuous luminosity in- 
jection, and the second term takes into account radiative 
energy losses in the blast wave. For q 
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t > to (Cohen & Piran 1999). Here to is a characteristic 
timescale for the formation of a self-similar solution, which 
is roughly equal to the time for the external shock to start 
to decelerate, and Ei mp is a constant which describes the 
impulsive energy input when q' > — 1 — k'. This is ob- 
tained under the assumption that a self-similar solution 
exists at t > to, hence it cannot be extrapolated down to 
t = 0. For t > to the bulk Lorentz factor of the fireball 
scales with time as T 2 oc t~ m , with m and n' connected 
by k! = m — 3 (Cohen, Piran & Sari 1998), and m — 3 for 
the adiabatic case (Blandford & McKee 1976). 

In the observer frame the time T is related to the fixed 
frame t (for which dr = cdt) by dT = (1 - (3)dt ~ dt/2T 2 , 

and T = f*(2T 2 )- 1 dt ~ t/[2(m + 1)T 2 } when t > t . The 
differential energy conservation relation can be now ex- 
pressed as dE/dT = L (T/T ) q - n(E/T), and the inte- 
grated relation is 



E = 



hj+i 



(C^ + ^mp^)^, T>T . (1) 

Here L = Lo(T/To) q is the intrinsic luminosity of the cen- 
tral engine, To = to/[2(m' + l)T 2 (to)] where m! is the self- 
similar index for t < to (usually ml = for the coast- 
ing phase), L = 2T 2 (t )C , q — (q' — m)/(m + 1), and 
K = K'/(m+l). Since K + q+l = (k' + q' + l)/(m + l), the 
comparisons between q and — 1 — k in the following dis- 
cussions are equivalent to the comparisons between q' and 
-1 - k' (Cohen & Piran 1999). Setting T = T , the total 
energy at the beginning of the self-similar expansion is the 
sum of two terms, Eq = LqTq/(k + q + 1) + E lmp . The first 
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term, for q > — 1 — k, is the accumulated energy from the 
continuous injection (with radiative corrections) before the 
self-similar solution starts (note that for q < — 1 — n the 
two terms no longer have a clear physical meaning since 
the first one is negative). The second term, -Ei mp , is the 
energy injected impulsively by the initial event. 

At different times the total energy of the blast wave 
given by Eq.(|lJ) may be dominated either by the continu- 
ous injection term (oc T' 9+1 )), or by the initial impulsive 
term (oc T~ K ). Which of these two is dominant at a par- 
ticular observation time T depends both on the relative 
values of the two indices (q + 1 and — k, Cohen & Piran 
1999), and also on the values of Lq and £i mp (Dai & Lu 
1998a, b). We then have three regimes: i) If q < — 1— k, the 
second term in ([!]) always dominates since the first term is 
negative. The fireball is then completely analogous to the 
impulsive injection case, ii) If q = — 1 — k, the solution ([!]) 
is no longer valid, and there is no self-similar solution, iii) 
If q > — 1 — k, the first term in ([lj) will eventually domi- 
nate over the second term after a critical time T c , and it 
is this term that will exert a noticeable influence on the 
GRB afterglow light curves. This latter case is of the most 
interest here. 

For a fireball blast wave decelerated by a homogeneous 
external medium with particle number density n, the en- 
ergy conservation equation at time t = r/c is 
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r 3 nm p c 2 T 2 = ^-{ctfnm p c 2 T 2 , t > t , (2) 



where r is the radius of the blast wave, and all other sym- 
bols have their usual meanings. This relation holds also 
for more general cases with a non-constant E. For the 
continuous injection dominated case, the energy E in (J2J) 
should have the same time dependence as the first term in 
the RHS of Eq. (IT]), giving T^ 1 oc t 3 T 2 , or 
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(3) 



Since in general T oc t~ m ' 2 oc r~ m l 2 oc T~ m / 2 ( m+1 ) and 
r oc 7 ll /( m + 1 ) j using Eq. (^|) the dynamics of the continu- 
ous injection-dominated case is 



T oc r"( 2_<? )/ 2 ( 2+ 9) oc jM 2 -^/ 8 r oc T (2+q ^ 4 



(4) 



For such a continuously-fed fireball, a forward shock propa- 
gating into the external medium and a reverse shock prop- 
agating back into the relativistic fireball will co-exist on ei- 
ther side of the contact discontinuity. The latter may per- 
sist as long as a significant level of energy injection is going 
on. The radiation spectra from these shocks are compli- 
cated by many factors. Here we address only the sim- 
plest case of the standard adiabatic external shock after- 
glow scenario, and assume that the reverse shock is mildly 
relativistic, as in the refreshed shock scenario. Follow- 
ing Meszaros & Rees (1997), Kumar & Piran (2000) and 
Sari & Meszaros (2000), one can work out the relationship 
between the temporal index a and the spectral index (5, 
where F v oc T a v$ . For the forward shock, the synchrotron 
peak-frequency i/£ oc YB'^ 2 m oc T 4 oc r 2m /( m+1 ) oc 
T -( 2 -g)/2 5 and the peak flux F f^ K (fY 5 ){n' e B'r/T) oc 

T 3p8 K T (3-m)/(l+m) x yl+g^ g^at a f = {2m(3 i + 3 - 

m)/(l + m) = (1 — q/2)[3^ + 1 + q. For the reverse shock, 



v T m = v f m /T 2 oc T 2 oc oc r-(2-3)/4 and F r _ 

TFf oc T 3 r 9 oc T( 6 - 3m )/ 2 ( 1+m ) oc T^ 9 ^/ 8 , so that 
a r = (2r7i/3 r + 6-3m)/2(l + m) = [(4 - 2q)fV + 9q + 6]/8. 
At any time the emission at a given frequency, may be 
dominated either by the forward or the reverse shock. The 
above scalings are derived for the slow cooling regime (Sari, 
Piran & Narayan 1998), which is usually satisfied. The fast 
cooling will change the scaling laws, but does not change 
the qualitative "index switching" picture proposed in this 
paper. All the above scalings reduce to the standard adi- 
abatic case by setting m = 3 and q = — 1. 

The injection-dominated regime begins at a critical time 
T c defined by equating the injection and energy loss terms 
in®, 



T r = Max^ 1, 



(K + q+1) 



Mrap 

LqTq 



To, (5) 



where T c > To ensures that a self-similar solution has 
already formed when the continuous injection law dom- 
inates. If initially the continuous injection term is more 
important, i.e., LqTq > E{ mp , then T c ~ To, and the dy- 
namics is determined by the continuous injection law as 
soon as the self-similar profile forms. However, if initially 
the impulsive term dominates (LqTo <SC -Eimp), the criti- 
cal time T c after which the continuous injection becomes 
dominant could be much longer than To, depending on the 
ratio of E- lmp and LqTq. 

Central engines with a continuous injection may, in addi- 
tion, have another characteristic timescale T, e.g. at which 
the continuous injection power law index (say q\ > — 1 — «) 
switches to a lower value q^ < — 1 — «. It is only for T > T c 
that the continuous injection has a noticeable effect on the 
afterglow light curve. Different central engines may have 
different values of T. In the following, we investigate the 
conditions under which a continuous injection can influ- 
ence the dynamics of the blast wave, and we consider the 
specific case of a millisecond pulsar as a central engine with 
a continuous injection following after the initial impulsive 
phase. 

3. HIGHLY MAGNETIZED MILLISECOND PULSAR AS THE 
CENTRAL ENGINE 

Although a wide range of GRB progenitors lead to a 
black hole - debris torus system (Narayan, Paczyhski & Pi- 
ran 1992; Woosley 1993; Paczyhski 1998; Meszaros, Rees 
& Wijers 1999; Fryer, Woosley & Hartmann 1999), some 
progenitors may lead to a highly-magnetized rapidly ro- 
tating pulsar (e.g. Usov 1992; Duncan & Thompson 1992; 
Thompson 1994; Yi & Blackman 1998; Blackman & Yi 
1998; Kluzniak & Ruderman 1998; Nakamura 1998; Spruit 
1999; Wheeler et al. 2000; Ruderman, Tao & Kluzniak 
2000). During the early stages of the evolution of these 
objects, the luminosity decay law could be very compli- 
cated. On the longer (afterglow) timescales that we are 
interested in, some short-term processes, such as the de- 
cay of the differential-rotation-induced toroidal magnetic 
field energy (Kluzniak & Ruderman 1998; Ruderman et al. 
2000), are no longer important, and the energy injection 
into the fireball may be mainly through electromagnetic 
dipolar emission. The spindown of the pulsar may also be 
influenced by gravitational radiation. 
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Assuming that the spindown is mainly due to elec- 
tromagnetic (EM) dipolar radiation and to gravitational 
wave (GW) radiation, the spindown law is — 7f2f2 = 
(B 2 .R 6 f} 4 )/(6c 3 ) + (32G/ 2 e 2 f2 6 )/(5c 5 ) (Shapiro & Teukol- 

sky 1983), where il and fl are the angular frequency and 
its time derivative, B p is the dipolar field strength at the 
poles, I is moment of inertia, R is stellar radius, e is el- 
lipticity of the neutron star. This equation can be solved 
for U as a function of T, with initial conditions Q = Do, 
Si = S7o f° r T — 0. The decay solution f2(T) includes both 
EM and GW losses, but the corresponding energy input 
into the fireball will be due to the EM dipolar emission 
only, i.e. L(T) = [B$R 6 Sl(T) 4 ]/(6c 3 ), which usually will 
not be a simple power-law. However, at different times 
the spindown will be dominated by one or the other loss 
terms, and one can get approximate solutions. When EM 
dipolar radiation losses dominate the spindown, we have 
O = Qo(l + T /Tem)^ 1 / 2 , or approximately £1 = £lo for 
T < T em , and 0. = ft (T/T em )~ 1/2 for T > T em . Here 



T — 



J^Lg = 2.05 x 10 3 s h 5 B~lP 2 _ 3 R-' 



(6) 



is the characteristic time scale for dipolar spindown, 
-Bp, 15 = £? P /(10 15 G), and Po,-3 is the initial ro- 
tation period in milliseconds. When GW radiation 
losses dominate the spindown, the evolution is Q ~ 
n (l + T/T gw )- 1 / A where T gw = (5c 5 /128G/e 2 ^) = 
0.91 s 2^ Pq _ 3 (e/0.1) -2 . GW spindown is important only 
when the neutron star is born with an initial f2o,>r2* ~ 
10 4 s" :L (e.g. Usov 1992; Blackman & Yi 1998). In such 
cases the ellipticity is large (e ~ 0.1) due to rotational in- 
stability, and the timescale for the GW-dominatcd regime 
is short, so that will be damped to below fi* promptly 
in a time % = [(fi /^*) 4 - tfTgw After Cl < fi,, GW 
losses decrease sharply, and the spindown becomes domi- 
nated by the EM losses. If the neutron star is born with 
flo < fi*, the spindown will be always in the EM regime, 
since the typical spindown time for GW radiation is much 
longer. 

For these regimes, the continuous injection is: 
(A) SIq < EM loss dominated regime. In this case 



where T em is given by (|6|), and 




T<T CT 
T > 71, 



(7) 



J em.Q 



- 1 c n 



1.0 x 10 49 erg S - x Bl 15 P^ 3 Rl (8) 



(B) f2 > initially GW dominated regime. The 

continuous injection luminosity can be divided into two 
phases, i.e., L = L em o/(l + T/T gw ) for T < % or 
L = L Pm .*/[l + (T-%)/T em ^] 2 for T > %, where 

m 2 /2T em ., ~ 1.0 x 
* = 3c 3 //P 2 i? 6 f> 2 ~ 



L emfi is given by (§), and L 



10 erg s 



1 d2 p 



3 i?g, where T el 



2.05 x 10 3 s I 4B B-=f 5 J? i _ 3 V. 

The above continuous injection luminosities have, for 
certain times, a temporal index g = > — 1, which may 



dominate the blast wave dynamics. The typical duration 
times for this flat injection law are T em for case (A) or T gw 
and T emt * for case (B). In case (B) there are two luminosity 
"plateaus" . The former is usually much shorter than T c 
unless a very dense medium is assumed (see below) , so it is 
unlikely to detect such a "two-step" injection-dominated 
case. If is close to fio, the second timescale for case 



(B), T eT 



may not be much different from T Prn . In the 



following, for simplicity, we discuss case (A) only, keeping 
in mind the possible extra complexity which case (B) may 
introduce in some extreme cases. 

During the time interval T c < T < T em one can expect 
a distinctive pulsar feature to show up in the lightcurve. 
We take as an example the behavior (0), setting q = 
in this regime, and q = — 1 otherwise (since the second 
slope q = — 2 would mimic in effect the standard impulsive 
q = — 1, m = 3 adiabatic case). For the slow-cooling exter- 
nal shock scenario (§2), the temporal decay index changes 
at T c from ot\ to ct2 = (2/3)ai + lj and changes back to 
ai after T em . The temporal index is related to the spec- 
tral slope through a,\ = (3/2)/? for the forward shock, 
and a\ = (6/? — 3)/8 for the reverse shock. This implies an 
achromatic bump in the light curve, which could provide a 
signature for a pulsar. However, not all millisecond pulsars 
would give such observable signatures. The condition for 
detecting this feature is T em > T c , which constrains the 
pulsar initial parameter phase space. Let us specify the 
pulsar case in Eq.(^), i.e., q = 0, and further assume K = 0, 
so that Eq.(|) is simplified to T c = Max(l, E lm p/LoTo)T . 
This gives two possibilities: 

(I) Eo is mainly due to the continuous injection, i.e., 
Ei m p <^ LqTq. This is the case considered in many 
pulsar GRB central engine models (Usov 1992; Dun- 
can & Thompson 1992; Thompson 1994; Spruit 1999). 
We then have T c = Tq and Eq ~ L em ^To (a fac- 
tor of 2 if the continuous injection and the prompt 
injection energies are comparable). Solving Eq = 
(47r/3)(2cr§T ) 3 nTOpC 2 r§ (in which m' = has been 
adopted), where T = T(T ) = E /AMc 2 is the ini- 
tial bulk Lorentz factor of the blast wave, we have T c = 
T ~ 0.33s Bp,i 5 P :!3i?^(ro/300)- 4 n- 1 / 2 . The condition 
T em > Tq then implies 



B 



p,15 



< 18AR 



.4/3 



T 1/3 R- 



(r /300) 4/3 n 1/6 . (9) 



(II) Eq is dominated by the impulsive term E- lmv . This 
could be the case if the central engine is a pulsar, the 
initial impulsive GRB fireball being due, e.g. to vv an- 
nihilation (e.g. Eichler et al 1989), dissipation of initial 
differential rotation (Ruderman et al. 2000), phase con- 
version from a neutron star to a strange star (Dai & Lu 
1998b), etc. In this case, since Ei mp > LqTo, we have 
T c = E imp /L ,em = (2E imp /IQl)T em . The condition 
T em > T c is simply E imp < (l/2)JOg, or 



Po,_ 3 < lAll^E.21%. 



(10) 



Since both T em and T c are large in this case, the contin- 
uous injection term may dominate at a later time. An 
additional constraint, to avoid rotational break-up of the 
pulsar, is 

P > Po(min) . (11) 
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The equations (H), ( |10[ ) and ([11]) define a region of the 
pulsar Pq,B p q initial parameter space (see Fig. 1), in- 
side which the continuous injection has the observational 
signature discussed above, i.e. a bump in the light 
curve. The separation between the two regimes (I, II) is 
defined by the condition L em-0 T ~ lO 51 ^, 



a P) 5i erg, or 



B, 



p,15 



6.7P 2 



-3^6 



^To/SOO^^n 1 /^-^ . Above this 



line (regime I), T c ~ T and one expects only one change, 
starting with the flat regime and changing to steep; be- 
low this line (regime II) one expects two changes, starting 
with the steep decay, followed by a flat regime, and a final 
resumption of the steep decay. 
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Po diagram for the initial parameters of a pul- 
sar born in a GRB. The enclosed areas (regime I and II) are the 
phase spaces where an achromatic pulsar signature is expected in 
the GRB afterglow lightcurves. Dotted-dashed lines denote various 



T ern expected. E- ln 



Po( mm ) = 0.5 ms have been adopted. 



10 bi erg s" 1 , To = 300, n 
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4. DISCUSSION 

If the central engine of GRB are fast-rotating pulsars, 
the afterglow lightcurves may show a distinctive, achro- 
matic feature, for pulsars whose initial parameters are 
within "pulsar signature region" region defined by equa- 
tions (JoJ) , (|ic| ) and (|lll). In this case the afterglow light 
curves flatten after a critical time T c (eq.§), and steepen 



again after a time T em (eq.||). This region of pulsar pa- 
rameter phase space includes "magnetars" or ultra-high 
field pulsars, which play a large role in some GRB pulsar 
progenitor models. 

The current data on early afterglows are insufficient to 
provide good tests for this feature. However, an inter- 
esting possibility is the achromatic bump observed in the 
afterglow of GRB 000301C, which for a limited time de- 
viates significantly from the standard broken power-law 
fit (Masetti et al. 2000). Possible explanations include 
running into a non-uniform ambient density (Berger et 
al. 2000), and modification by a microlensing event (Gar- 
navich, Loeb & Stanek 2000). We suggest here a third 
possibility, that the bump may be caused by the pulsar 
signature discussed above. Taking (Berger et al 2000) 
ot\ ~ —1.28 for the principal temporal index (before the 
bump and before the decay ascribed to a jet transition), 
the temporal index during the pulsar signature is expected 
to be «2 ~ 0.15, which seems reasonable to fit the achro- 
matic bump. Since the feature occurs a couple of days 
after the GRB trigger, the initial pulsar parameters are 
in regime II. Taking -Ei mp .5i ~ 1.1 from the observations 
(Berger et al. 2000), the pulsar parameters follow from 
setting T em ~ 3.8 d and T c ~ 2.5 d, which gives Pq ~ 3.4 
ms and B p ^ ~ 0.27 (Fig.l). Detailed x 2 fit to the data 
would be necessary to validate this proposal. 

Another relevant observation may be the recent Fe line 
detection in the X-ray afterglow of GRB 991216 (Piro 
et al 2000), which in one interpretation would require a 
continuously-injecting central engine (Rees & Meszaros 
2000). If L em ~ 10 47 erg s _1 is assumed 37 hours after 
the burst (Rees & Meszaros 2000), this continued injection 
could be due to a pulsar with B p ^ 5 ~ 0.15 and P < 1.2 
ms, which would imply a signature bump at ^1 hour. 
Such a bump is not seen in this afterglow (e.g. Halpern 
et al. 2000). However, a slightly weaker luminosity (e.g. 
Lem < 3 x 10 46 erg s -1 ) could also explain the Fe line fea- 
tures, by assuming a slightly larger Fe abundance (which is 
very low in this model) . The required pulsar can therefore 
be more magnetized, and the characteristic times for the 
signature bump would be expected early enough to have 
evaded detection in this GRB. 

In conclusion, the detection of or upper limits on such 
characteristic afterglow bumps by missions such as HETE2 
or Swift may be able to provide interesting constraints on 
magnetar GRB models and their progenitors. 
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